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Through detailed radiative transfer modeling, we present a disk+cavity model 
to simultaneously explain both the spectral energy distribution (SED) and Sub- 
aru if-band polarized light imaging for the pre-transitional protoplanetary disk 
PDS 70. In particular, we are able to match not only the radial dependence, but 
also the absolute scale, of the surface brightness of the scattered light. Our disk 
model has a cavity 65 AU in radius, which is heavily depleted of sub-micron- 
sized dust grains, and a small residual inner disk that produces a weak but still 
optically thick Near-IR excess in the SED. To explain the contrast of the cavity's 
edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized 
dust inside the cavity is required. The total dust mass of the disk may be on 
the order of 1O _4 M , only weakly constrained due to the lack of long wavelength 
observations and the uncertainties in the dust model. The scale height of the 
sub-micron-sized dust is ~ 6 AU at the cavity edge, and the cavity wall is opti- 
cally thick in the vertical direction at if -band. PDS 70 is not a member of the 
class of (pre-)transitional disks identified by Dong et al., whose members only 
show evidence of the cavity in the millimeter-sized dust but not the sub-micron- 
sized dust in resolved images. The two classes of (pre-)transitional disks may 
form through different mechanisms, or they may simply be at different evolution 
stages in the disk clearing process. 

Subject headings: protoplanetary disks — stars: pre-main sequence — stars: 
individual (PDS 70) — radiative transfer — circumstellar matter 



1. Introduction 



Recent years have witnessed a rising interest in a special kind of protoplanetary disks, 
in which a giant central cavity is present and reveals itself in the spectral energy distribution 



(SED, e.g. 


Strom et al. 


1989 Skrutskie et al 


. 1990 Calvet et al. 


2005; 


Furlan et al. 


2006 


Espaillat et al. 2007a|b 


2008)), or in sub-mm interferometry observations (e.g. 


Pietu et 


al. 2006[ Hughes et al. 


2007 Brown et al. 


20091 Hughes et al. 


2009 


Isella et al. 


2010 


Andrews et al. 2011 Mathews et al. 2012; C 


Jieza et al. |2012b; Isella et al. 2012 


Andrews 


et al. 2012 


. Depending on whether a small residual optically thick inner disk is left at the 



center, these objects could be classified into two categories — pre-transitional disks (with an 



1 Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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optically thick residual inner disk, Espaillat et al. 2010), and transitional disks (without an 



optically thick residual inner disk, Espaillat et al.||2007b ). Studying these objects at multiple 
wavelengths is of great interests, because the inner disk clearing process may be signposts 



of planet formation and disk evolution in general ( 


Zhu et al. 


2011 


Williams & Cieza 


2011 


Dodson- Robinson & Salyk 2011 Birnstiel et al. 


2012 


Cieza et al. 


2012a 


)• 



So far, the study of (pre-)transitional disks was done primarily via measurements of the 
SED and resolved images at sub-mm wavelengths. More recently, direct imaging at optical 
to near infrared (NIR) started playing a crucial role. Protoplanetary disks usually have 
dust grains with sizes ranging from sub- micron (called small dust in this study), similar 
to the pristine interstellar medium (ISM) dust, to mm or larger (big dust), forming as a 



result of grain growth and coagulation (Dullemond & Dominik 2005 Birnstiel et al. 2012). 



The SED of a disk is degenerate with many disk parameters, since it is spatially integrated 
over the entire disk. Resolved images at sub-mm, on the other hand, provide more detailed 
information, but they are only sensitive to the distribution of big dust due to its dominance 
in opacity at these wavelengths. Using 8-m class ground based telescopes equipped with 
Adaptive Optics and/or coronagraph system, direct imaging at optical to NIR wavelengths 
provides resolved disk maps with high spatial resolution (~ 0.06") and small inner working 
angle (< 0.2"). These maps are directly related to the distribution of small dust. 

An ongoing survey scale project, the Strategic Explorations of Exoplanets and Disks 



with Subaru (SEEDS, Tamura 2009), is directly imaging a large sample of protoplanetary 
disks at NIR (J, H, and K) bands in a systematic way, using the High-Contrast Corono- 
graphic Imager for Adaptive Optics (HiCIAO, Suzuki et al.|2010 ). SEEDS disk observations 
specialize in taking polarized intensity (PI=a/Q 2 + U 2 , which Q and U are components in 
the Stokes vector) images of disks, which greatly enhances our ability to probe disk struc- 
ture (especially at the inner part), by utilizing the fact that the central source is usually not 
polarized, so that the stellar residual in PI images is much smaller than in full intensity (FI) 
images dPerrin et al.]2004{ [Himdey et al.|2009||Quanz et al.|2011| ). The signal in NIR imaging 
usually comes from scattering of starlight off the surface of the disk, since dust at separations 
of ~ 10 AU in the disk is not hot enough to contribute significantly at NIR wavelengths. 



A number of protoplanetary disks have been studied in SEEDS (LkCa 15, Thalmann et al. 



2010, AB Aur, Hashimoto et al. 2011 SAO 206462, Muto et al. 2012, MWC 480, Kusakabe 



et al. 2012, and UX Tau A, Tanii et al. 2012), and more are on the way. 



Recently, Dong et al. (2012) pointed out that in a large sample of (pre-)transitional 
disks whose central cavity have been confirmed in resolved images at sub-mm, SEEDS did 
not find the cavities in their NIR images. This discrepancy has been interpreted as evidence 
for the cavity existence only in big dust, with the small dust having a continuous distribution 
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in surface density from the outer disk all the way to the inner working angle of SEEDS 
(~ 071 — Of! 15, or about ~ 15 — 20 AU in Taurus). This may be explained by the dust 



filtration model proposed by Paardekooper & Mellema (2006); Rice et al. (2006); Zhu et al. 



(2012). In their model, the pressure bump in the disk acts like a filter, which filters through 



the small dust but traps the big dust, resulting in the depletion of big dust inside the cavity. 
In this article, we study the structure of the pre-transitional disk PDS 70 using detailed 



radiative transfer modeling, following the observational paper by Hashimoto et al. (2012). 



PDS 70 is a K5 type weak-lined T Tauri star located at the Centaurus star-forming region 
(coordinate 14 08 10.15 -41 23 52.5), with an age < lOMyr and a distance estimated to be 
~ 140pc ( |Gregorio-Hetem fc Hetem||2002| |Riaud et al. | |2006[ |Metchev et aL|2004t |Hashimoto 
et al. 2012). The photometry of this star has been measured at multiple wavelengths (e.g. 



Metchev et al. 2004 Kessler-Silacci et al. 2005), and its disk has been detected in scattered 



light by VLT if-band imaging (Riaud et al. 2006). SEEDS directly imaged PDS 70 in H 



band using the polarization differential imaging (PDI) mode on Feb. 27, 2012. The Subaru 



observation and the data reduction were described in detail in Hashimoto et al. (2012, which 



also presented the Gemini L'-band imaging). Unlike most previous objects of this type, the 



cavity in PDS 70 is first detected by SEEDS NIR imaging. And unlike the objects in Dong 



et al. (2012), it indeed has a cavity in the small dust, which also contradicts the predictions 



in various dust filtration models mentioned above. 

Radiative transfer modeling of (pre-)transitional disks has been mostly used to fit the 



SED and sub-mm observations in the past (e.g. Hughes et al. 2009 Andrews et al. 2011 



Mathews et al. 2012 Andrews et al. 2012). Here we demonstrate the power of modeling 



in NIR, by producing synthetic disk images in scattered light and comparing them with 
observations. Sensitive to dust at different sizes and locations, SED, sub-mm observations, 
and NIR imaging can be used to probe the disk structure in different and complementary 
ways. Particularly, NIR imaging can provide effective constraints on the scale height of the 
disk, and the depletion factor for the small dust inside the cavity, neither of which is well 



constrained using the other two types of observations (Andrews et al. 2011; Dong et al. 



2012). 



This paper is the first in a series, in which we study the structures of (pre-)transitional 
disks by modeling SEEDS NIR imaging data and observations at other wavelengths. The 
paper is arranged in the following way. In Section [2] we introduce the method for our 
radiative transfer modeling. The fiducial model of PDS 70 is presented in Section [3j where 
we simultaneously fit both its SED and scattered light image. The constraints on various 
disk+cavity parameters from current observations are discussed in Section |4j where we focus 
on the ones which are directly probed by NIR imaging and are important in revealing the 
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disk evolution. We discuss the possible formation mechanisms of (pre-)transitional disks in 
Section |5j followed by a short summary in Section [6] 



2. Radiative transfer modeling 

We use the Monte Carlo radiative transfer (MCRT) code developed by B. Whitney 
(2012, in prep, see also |Whitney et al.|2003a|b) |Robitaille et~aL|2006[ ) to model the SED and 
SEEDS H-b&nd PI imaging of PDS 70. The SED data collected from literature, the SEEDS 



observation, and its data reduction are described in detail in Hashimoto et al. (2012). MCRT 



simulations are run with 4 x 10 7 photons for high S/N images. Disk setup is largely adopted 



from the model in Dong et al. (2012), which is briefly summarized here. We construct an 
axisymmetric disk 200 AU in radius on a 600 x 200 grid in spherical coordinates (R, 8), where 
R is the radial component and 9 is the poloidal component (6 = 0° is the disk mid-plane). 
We uniformly deplete the central region of the disk to form a cavity with radius -R cav , i-e. 
reduce dust surface density by a constant factor within this radius. The disk is assumed to be 
at ~140 pc. The inner radius of the disk (within the cavity) is self-consistently determined 
by dust sublimation (T sub ~ 1600 K). The accretion rate of PDS 70 is considered to be very 
low, since the object has been identified as a weak line T Tauri object with an Ha equivalent 
width of 2A ( |Gregorio-Hetem fc Hetem||2002[ ). For the central source, we assume a pre-main 
sequence star of spectral type K5, radius 1.35 R & , mass 1.0 M , and temperature 4500 K, 



as suggested by Gregorio-Hetem & Hetem (2002). The surface density in the outer 

disk is taken as 

E(R) = E ^e- R / R % atR> i? cav , (1) 



where R c is a characteristic scaling length assumed to be 50 AU, and the normalization So is 
determined by the total dust mass of the disk M dust . The choices of the surface density radial 



profile and R c are common in protoplanetary disk modeling (Williams & Cieza 2011 Andrews 



et al.|20iT ), and represent a snapshot of a solution for a fully viscous disk with a constant a 
parameter (as in Shakura & Sunyaev 1973 a disk model) and T oc i? -1 / 2 (Hartmann et al. 



1998). Inside the cavity, the disk is uniformly depleted, with surface density going as 



R 



-R/Rc 



, at R < R c 



(2) 



where 5 cav is the constant depletion factor (which may be different for different dust popu- 
lations). 

Various grain evolution models predict that the pristine dust grains in the disk coagulate, 



grow to bigger size, and subsequently settle to the disk mid-plane (Dullemond & Dominik 
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2004a[b| |2005[ ). To take this effect into account, we assume a two component model of dust 
distribution: a thick disk with small grains (sub-micron-sized), and a thin disk with large 



grains (up to ~mm-sized). Below we will use subscripts "b" a 
relating to big and small dust, respectively, and superscripts 



for various quantities 
'c" to indicate the 



outer disk and the cavity, respectively (for example, T, c s represents the surface density of the 
small dust inside the cavity). Both big and small grains are assumed to have a Gaussian 
density profile in the vertical direction, 



p(R,z) 



E(R) 



2 /2h 2 



27ch 



(3) 



where h is the scale height, with h\> (scale height of the big dust) assumed to be much 
smaller than h s (scale height of the small dust). We note that the scale heights are provided 
to the code as input parameters, instead of being self-consistently determined from the disk 
temperature calculated in the simulations. In our fiducial model shown below, we check this 
assumption and find that the input is consistent with the output (Section [3]). Radially, both 
scale heights vary with radius as 

h oc R^ (4) 

where /3 is the constant power law index. Our experiments show that as long as the big 
dust is settled to the disk mid-plane (i.e. h h <C h s ), the details of its vertical distribution 
hardly affect the details of NIR images and SED, while both crucially depend on the spatial 
distribution of the small dust, as we will show below in Section [4~3 



The total mass of the big dust in the disk is / x Md us t, where / is a variable parameter. 
For the small dust, we assume a size distribution as in the standard ISM dust model from 



Kim et al. (1994) up to a maximum size of 0.2 /zm (roughly a power-law size distribution 
n(s) oc s p with power index p ~ 3.5). The dust composition is assumed to be 70% mass in 



silicate and 30% mass in graphite, with the properties for both adopted from Laor Sz Draine 
(1993). Preferentially forward scattering (i.e. Mie scattering) is assumed. For the big dust 



we primarily use Model 2 from Wood et al. ( 2002[ ) (from now on called steep-big-dust), which 
assumes a power-law size distribution n(s) oc s p with p = 3.5 up to a maximum size of 1 mm, 



and a composition of amorphous carbon and astronomical silicates. In Section 4.2 Model 3 



from Wood et al. (2002) is also tried (from now on called flat-big-dust), whose only difference 
from the steep-big-dust model is that it assumes p = 3.0 ("flatter" size distribution, so more 
mass at the large size end). As a consequence, the flat-big-dust model has lower opacity 
than the steep-big-dust model at wavelengths A < 1mm. The opacity of these dust models 
are shown in Figure [TJ 

PDS 70 shows a cavity with radius ~ 60 — 70AU, and the disk is inclined by about 



45° — 50° (Hashimoto et al. 2012). In our fiducial model shown below we choose 65 AU 
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0?46) and 45° for the two, as they produce model SED and image which match observations 
reasonably well (Section [3]). We note that a mild deviation from these "fiducial" values can 
be tolerated without much difficulty (i.e. ±5 AU in cavity size and ±5° in geometry). These 
parameters are fixed below to keep the models simple. 

We produce if-band polarized images from MCRT simulations. To obtain realistic 
images which can be directly compared with SEEDS observations, the raw simulated images 



are post-processed, as described in detail in Dong et al. (2012). The SEEDS observations 



of PDS 70 were conducted in PDI mode without a coronagraph (Hashimoto et al. 2012) 



so that we produce mock images in the same mode in this work. The raw model images 
of the entire disk+star is convolved by a real if-band SEEDS point spread function (PSF). 
The inner working angle of the images, ip- m , is assumed to be / .'15 in radius (~21 AU at the 
estimated distance of PDS 70). We measure the surface brightness radial profile (SBRP) 
of the disk along major and minor axes, by calculating the average SB of the pixels within 
±22.5° on each side of the axes at various radial bins (T/05 in width, the same as we measured 



for the SBRP of the SEEDS PDS 70 image (Hashimoto et al. 2012). 



3. The fiducial disk model for PDS 70 



Here we present a fiducial disk model for PDS 70 to fit all the observations. The disk 
parameters in this model are listed in Table [TJ and the surface density of the small dust is 
plotted in Figure [2j In summary, the disk has a heavily depleted cavity whose radius is 65 
AU i? cav , with a depletion factor of 1000 for both the big and small dust (5 cav ,b and £ C av,s? we 
note that only <5 caVjS is constrained by the current observations, not £ cav ,b, see below). The 
inner-most disk (on AU scales) and the cavity wall are both optically thick in the vertical 
direction at 1 /xm, representing the peak of the stellar spectra. The vertical optical depth 
due to the small dust is ~ 14 at the inner edge of the disk (~ 0.06 AU), and ~ 4 at the cavity 
wall. On the other hand, the gap in between the inner rim and the outer disk (i.e. from ~1 
AU to .Rcav) is optically thin to stellar radiation, which justifies the classification of PDS 70 
as a pre-transitional disk (Espaillat et al. 2010). The total dust mass Md ust is 3 x 1O _5 M , 
and most of it is in the big dust, with a big-to-small-dust ratio at about 30:1 (i.e. / = 0.97). 



However, as we will discuss in Section 4.2 , the constraints on the big dust are rather weak 
(including the total amount, its structure, and the big-to-small-dust ratio), given the current 
data. We assume scale height power law index 1.2 for both dust populations (/3b and /3 S ), 
both inside and outside the cavity (indicated by subscripts "c" and "o"), reasonably close to 
the canonical value ~ 1.25 — 1.3 for irradiated disks ( Chiang fc Goldreich||1997 Hartmann et 



al. 1998). The gas in the disk is in hydro-static equilibrium in the vertical direction, so that 
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/i gas is set by the vertical temperature profile, through h gas « c s /f2 and c s ~ ^Jk B T / Ji, where 
c s is the sound speed, Q is the rotational angular velocity, T is the mid-plane temperature, 
and \x is the average molecular weight. Various grain settling models predict that the small 



grains tend to be well coupled with the gas (Dullemond & Dominik 2004b), as a result the 



two would share similar vertical distribution, h s ~ /i gas . The small dust scale height is 6.0 
AU at R cav , both inside the cavity (h^) and outside the cavity (h°), which corresponds to a 
mid-plane temperature of 28. 2K. This value is close to the output from the radiative transfer 
calculation, 31. 8K. The big dust scale height inside and outside the cavity (h^ and h^) is 1/5 
of the value for the small dust. 

The model H band PI images are shown in Figure [3j along with the observed SEEDS 
image. The SED and SBRP (defined in Section [2]) of the image are shown in Figure |4j 
Our model SED matches the observations very well. For the scattered light image, our 
fiducial model looks very similar to SEEDS image, both revealing a clear cavity with size 
~ 0'. / 5 on the major axes. If fit by an ellipse, the disk center has an offset ~ 9 AU from 
the star, roughly along the minor axis, and towards the far side of the disk, which is due to 
the back illumination of the wall. This offset is consistent with observation (the measured 



offset in SEEDS image is ~ 6AU in the same direction, see Hashimoto et al. 2012). The 
small (bright) structures that appear close to the inner working angle in the SEEDS image 
is probably artificial, mostly likely caused by observational noise or the stellar residual in 
polarized light. The right side of the Subaru image is slightly brighter (in scattered light) 
than its left counterpart, which is not reproduced by our axisymmetric fiducial model. This 
asymmetry may be intrinsic, caused by differences in grain properties, or disk structure (i.e. 
scale height or the amount of dust), on the two sides. We note that since we focus on 
studying the global scale structure, particularly the general properties of the cavity (i.e. its 
size and depletion factor), using axisymmetric models, we do not address the local details 
and non-axisymmetric structures in this work. 

Quantitatively, we reach good agreement with the measurement of the SBRP along the 
major axis. The only obvious deviation happens at large distance (>"1, well beyond the 
cavity edge), where unlike the model, the observed image flattens out to a (roughly constant) 
background noise. Due to the axisymmetry of our model, the SBRPs along both directions 
at its semi-major axis are the same, while observationally these profiles are slightly different, 
due to reasons discussed above. Here we emphasize that we achieve good agreement not 
only for the radial dependence (i.e. the slope), but also for the absolute scale of the surface 
brightness (the vertical axis in all the SBRP plots in this study is in actual physical units, 
and the curves were not rescaled). 

The agreement between our model and SEEDS observation on the minor axes is not 
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perfect, but nevertheless the two agree on major qualitative features: the surface brightness 
on the far side of the disk peaks at a larger radius, and it decreases outward slower than on the 
near side. The flux in polarized scattered light is determined by the product of polarization 
fraction, and the intensity of the full intensity (i.e. PI=(PI/FI) xFI). Along the major axis, 
the scattering angle is nearly 90°, which results in a maximum polarization fraction (PI/FI) 
due to the phase function of small dust. On the other hand, along the minor axis the 
scattering happens closer to the center, and at angles off 90°, resulting in small PI/FI. On 
the far side, both the upper and lower edge of the cavity are visible, so that the ring is wider; 
on the near side the lower edge is blocked by the outer disk, however forward scattering leads 
to a bigger FI there. We note that generally speaking, the comparison along the major axis 
is more valuable in constraining the disk+cavity structure, because disk is spatially more 
extended in this direction so that it is better resolved by observations with a fixed spatial 
resolution in all directions. Also, along the minor axis the features of interests (i.e. cavity 
edge) are present at closer separation from the center (also ipi Q ), where the photon noise is 
generally larger. 

We note that the parameters and geometry of our fiducial model come from an overall 
consideration of fitting the SED, NIR image, and radial profile of the scattered light simul- 



taneously, instead of simply measuring from the image as in Hashimoto et al. ([2012) . Also 



the fiducial model presented here is by no means unique, i.e. the only one which can provide 
a good fit to all the observations, as the constraints on some of the disk+gap parameters are 
rather weak (Section El). 



At last, we comment on the previous VLT i^-band imaging of PDS 70 reported by 



Riaud et al. (2006), and their derived disk model. Observationally, the inclination and 



position angle found by Riaud et al. (2006) is similar to our results (Hashimoto et al. 2012). 



However, our Subaru if-band PI images provide far more details and are at a higher quality 



than the VLT i^-band full intensity images presented by Riaud et al. (2006), which detected 



the disk in scattered light, but did not reveal a clear giant cavity structure, and showed a 



jet structure that is not present in our observations. Riaud et al. (2006) employed a smooth 



disk model with radius larger than 500 AU, and a total dust mass between 0.001 to 0.002 
M & to reproduce the observations, which is one order of magnitude lager than in our fiducial 
model. The difference is mostly caused by the fact that the inner region is heavily depleted 



in our models but not in Riaud et al. (2006). 
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4. Constraints on various disk and cavity parameters 

The fiducial model presented above provides a reasonably good fit to observations, 
and gives us a basic idea about the disk+cavity structure. In this section, we intend to 
determine the constraints on some of the parameters from modeling and fitting, so that we 
can understand what these observations are really telling us. In general, disk properties at 
short wavelengths (i.e. scattered light images and NIR excess in the SED) are sensitive to 
the spatial distribution of the small dust, while this dependence shifts to the big dust at long 
wavelengths (i.e. FIR excess and sub-mm observations), due to the difference in opacity of 
the two populations (Figure [T]). Since there are too many free parameters in our model, it is 
not realistic to vary every one of them and study their effects. Rather, we narrow our scope 
to few key parameters: 

1. The ones which are crucial in revealing the evolution of protoplanetary disk, and the 
formation mechanisms of (pre-)transitional disks. These are the depletion factor of 
both dust populations, and the total mass of the disk. 

2. The ones which can only be effectively constrained by scattered light images, such as 
the scale height and depletion factor of the small dust, which has a large degeneracy 
in SED and sub-mm observations. Since SEEDS virtually opens a new window to 
systematically study a large uniform sample of disks using scattered light images, we 
intend to provide an example to demonstrate the power of NIR imaging for probing 
disk structure. 



Below we study the role played by these factors both in disk SED and scattered light 
imaging, by exploring the corresponding parameter space around the fiducial model. We 



examine the cavity depletion factors for both dust populations (deb and des) in Section 4.1 
the total dust mass (Md us t) in Section 4.2 and the scale height of the small dust in the 



cavity and the outer disk [h c s , /3g, h°, and /3°) in Section 4.3 Models and their parameters 
in each section are listed and described in the corresponding block in Table [T] For the SED 
we focus on comparing with the photometry data, and for the scattered light we only look 
at the SBRP along the major axis. 



4.1. The depletion factor inside the cavity 

The cavity depletion factor of the small dust <5 caViS could be constrained by both the 
SED and the scattered light image. Qualitatively, depleting the cavity more in small dust 
reduces the short wavelength excess on the SED, and enhances the contrast of the cavity 
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in the image. Technically, the SED is more sensitive to the depletion in the inner part of 
the cavity (at AU scales) since most of the short wavelength excess is produced there, while 
the scattered light images are more sensitive to the depletion in the outer part of the cavity, 
where the contrast of the cavity edge is produced (the inner part of the disk within ip- m 
cannot be directly accessed in imaging observations). As a result, the two could in principle 
be constrained "independently". However, to simplify our discussion, we use a uniform 
instead of radius-dependent £ cav ,s ( an d <5 caVj b), an d note that this treatment does not affect 



our conclusion (see Dong et al. 2012 for a discussion on radius-dependent S c 



Figure [5^, shows the effect of 5 cav ,s on the SED (<5 C av,b is locked to <5 C av,s, but this hardly 
affects our result, as shown below). The surface density profile of the small dust (S s ) for 
these models is plotted in Figure [2j Whether 2 — 20/im excess strongly correlates with 
^cav,s sensitively depends on whether the innermost disk (at AU scale) is optically thick or 
not (characterized at the peak of the stellar spectra, ~ 1/xm). When the innermost disk is 
optically thick (pre-transitional disks), the IR excess is almost independent of the amount of 
small dust inside the cavity (the fiducial model and model SCM1). However, once S s at the 
inner disk decreases below the optically thick limit, and enters the transitional disk phase, 
the 2 — 20/im excess drops significantly as a result of the decreasing 5 C av,s (model SCM2 to 
SCM4). If the inner disk is completely depleted (SCM4), excess below ~10 fim disappears, 
and a nearly blackbody thermal component peaking at ~ 40/zm clearly reveals itself, which 
arises from the cavity wall. Though not as isolated as in SCM4, this wall emission signal is 
prominent in all models. Comparing with observation, we conclude that the NIR excess in 
PDS 70 is consistent with an optically thick innermost disk, which requires £ cav ,s > 10~ 4 . 

Figure [5)3 shows the SBRP of the convolved image for the above models. As expected, 
decreasing <5 caVjS makes the inner disk fainter. Moreover, the blocking effect due to the inner 
disk is reduced, so that a bigger area of the cavity wall is illuminated, and more starlight 
(rays closer to the disk mid-plane) reaches the wall and outer disk, increasing their brightness 
in scattered light. Similarly, increasing <5 caVjS tends to wipe out the signal of the cavity in 
scattered light (i.e. the bump at ~ 0?05), and makes the SBRP smoother. The observed 
SEEDS SBRP is broadly consistent with 5 caVjS ~ 10~ 3 , while a modest deviation (i.e. a factor 
of ~ 3) around this value could be tolerated without too much difficulty 

On the other hand, the cavity depletion factor of the big dust could only be directly 
constrained by sub-mm properties of the disk ( Andrews et al.|[20TT Cieza et aL]|2012b Isella 



et al.||2012 Mathews et al.||2012 ). Figures ^\p,d show that changing £ cav ,b hardly affects either 
the disk SED or NIR image (we fix the small dust to isolate the effect, i.e. the depletion 
factor for the small dust component inside of the cavity is fixed to <5 caVjS ~ 10~ 3 in all the 
models plotted here). We leave the constraint on 5 cav ,b to future sub-mm observations. 
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To conclude, 5 cav s ~ 10 3 is needed for the models to be consistent with both the SED 
and SEEDS image. The constraint is modest, roughly half a dec around the fiducial value. 



4.2. The total dust mass of the disk 

The total dust mass of the disk Md ust is a quantity crucial for a thorough understanding 
of disk evolution. In general, Mdust is determined by the big dust since it dominates in mass 
over the small dust, and for (pre-)transitional disks like PDS 70, most of the mass resides in 
the outer disk since the inner disk is heavily depleted. M^ ust is best constrained by sub-mm 
or mm observations, because the disk is usually optically thin at these wavelengths, so Md ust 
can be calculated from the measured sub-mm or mm flux, assuming a dust temperature 



and opacity model (i.e. Equation 2 and 3 in Williams & Cieza 2011). However, the total 



disk mass derived in this way normally contains large uncertainties, introduced by poorly 
unconstrained dust opacity and gas-to-dust-ratio, if converting dust mass to total gas mass 



(Panic et al. 2008). 



For PDS 70, we constrain the disk mass using the longest wavelength photometric data 
point available, which is at 160/aul. Our models are all vertically optically thin (sometimes 
only marginally) at 160/mi given our dust models. Models BOM1 and BOM2 (both have the 
steep-big-dust model) in Figure |6^i show the SED dependence on Md us t- While the fiducial 
choice of M dust = 3 x 1O~ 5 M agrees well with the observed photometry at 160/im, a factor 
of ~ 3 deviation in M^ ust from this value leads to a factor ~ 2 difference in 160/xm flux, 
while the SED at wavelengths shorter than ~ 100/im is largely unchanged. 

One issue which deserves special attention in this exercise of Md us t determination is the 
effect of the big dust model. Since the long wavelength flux in the optically thin regime is 
proportional to the opacity of the dust at that wavelength, changing the dust model for the 
big dust has a direct impact on M dust . Figure [6] shows two examples in which the flat-big-dust 
model is assumed instead of the steep-big-dust model. As mentioned in Section [2] and shown 
in Figure [TJ the flat-big-dust has a flatter size distribution, resulting in smaller opacity at 
wavelengths shorter than ~lmm. Model BOM3-flat has an identical set of disk parameters 
as the fiducial model. While the difference in SED between the two models is negligible at 
A < 100/im, the IR excess in BOM3-flat sharply drops below the fiducial model at longer 
wavelengths. To pull the 160/im flux back to the observed value, the disk mass needs to rise 
to 1.5 x 1O~ 4 M , as in model BOM4-flat. On the other hand, since all the BOM models 
have identical spatial distribution for the small dust, their scattered light images are almost 
the same, as shown in Figure [6b. 



-14- 



We conclude that the total dust mass of PDS 70 M dust is probably on the order of 
1O~ 4 M , depending on the dust model for the big dust: a flatter grain size distribution for 
the big dust corresponds to a larger Mdust- Observations at longer wavelengths are needed 
to distinguish different big dust models and to pin down the disk mass. 



4.3. The scale height of the small dust 

The scale height of the small dust h s , which dominates the absorption of starlight, is 
a central quantity involved in determining almost all the observable properties of the diskj^] 
Scale height and total mass of the small dust determine the shape of the disk surface, on 
which the stellar radiation is absorbed and re-emitted at longer wavelengths. Part of that 
emission escapes and gets observed, and part of it goes deeper into the disk interior to 
heat the disk, subsequently being reprocessed into disk emission at even longer wavelengths. 
Despite its importance, unfortunately very few types of observations can determine h B in a 



straightforward manner (see the example in|Muto et al. 2012, where the shape of the spiral 



waves in SAO 206462 is used to estimate h s ). However, scattered light is an excellent tool 
for probing h s , since its signal comes right from the disk scattering surface, which is directly 
determined by h s . 

In this section, we use several models to investigate the effects of h s on disk SED and 
NIR image. Figure [7^,b show the result of models SOS1-4, in which the scale height at the 
outer disk h° is varied. We find that the radial dependence of h° has a minimal effect on 
both the SED and scattered light image (models SOS1 and SOS2), while the overall scale 
of h° plays a prominent role in both the SED and NIR image. A much larger area of the 
cavity wall in model SOS4, whose /i° is 1.5 times the fiducial value, is directly exposed to 
the stellar radiation, due to its larger height. As a result, its SED is a factor of ~ 2 higher at 
the peak of the wall emission (~ 40//m) compared with the fiducial model and the observed 
SED, and the cavity edge is ~ 50% brighter. 

Figures [7^, d show the effect of h s inside the cavity (h^, model SCS1-4). Similar to h°, 
the absolute scale of the scale height plays a bigger role than its radial dependence. Model 
SCSI and SCS3, both having a thinner innermost disk, produce less NIR excess and slightly 
less scattered light from the inner disk. Also, since less starlight is blocked by the innermost 
disk and more of it reaches the cavity edge, the two have more MIR excess around 40/im 



2 The big dust is generally considered to have settled to the disk mid-plane, and its opacity at the peak of 
the stellar spectrum is much lower than that of the small dust. As a result, the detailed vertical distribution 
of the big dust does not have a prominent effect on the SED and the scattered light image. 
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and a slightly brighter cavity edge at //-band. On the contrary, model SCS4 has a thicker 
disk inside the cavity, resulting in more short wavelength excess and less long wavelength 
excess, due to less starlight reaching the outer disk. Its inner disk is also brighter and outer 
disk is fainter in the scattered light image. In general, the effect of on the SED is modest, 



echoing the finding in Section 4.1, that once the inner disk becomes optically thick, the NIR 



excess is insensitive to the distribution of small dust. 

In sum, the scale height of the small dust at R c &v, h/R ~ 0.09, is relatively well deter- 
mined by both the SED and the scattered light image (note that this value is consistent with 
the output temperature from radiative transfer calculation, Section [3]). On the other hand, 
the constraint on the radial dependence of h s is weaker, as a broad range of (5 from 1.15 to 
1.25 (both inside and outside the cavity) in the parameter space we explored here does not 
contradict the observations. 



5. Indication on the formation of (pre-)transitional disks 



PDS 70 is a special pre-transitional disk, in the sense that unlike most of its previous 
cousins, whose cavities were usually first inferred from the shape of the SED and then 
confirmed by resolved sub-mm images, the cavity in PDS 70 was first found in SEEDS 



scattered light image (see also Honda et al. 2010 for AB Aur). Follow up interferometer 
observations at long wavelengths, for example using the Submillimeter Array (SMA) or 
the Atacama Large Millimeter Array (ALMA), are needed to check the existence of the 
cavity in the big dust. If it is indeed confirmed, then this object will represent a class 



of (pre-)transitional disks which are in clear contrast to the ones discussed in Dong et al. 



(2012). Dong et aL](2012) reported the discovery of a class of (pre-)transitional disks, whose 



cavities are confirmed using the SMA interferometer but not seen in SEEDS scattered light 
images, despite the fact that the inner working angle of the SEEDS images is small enough 
to reveal these cavities (if they exist). In addition, observations at long wavelengths will 
help to determine the total dust mass of the disk, which is currently poorly contained to 
be ~ 1O~ 4 M due to the lack of long wavelength observations and the uncertainties in dust 
model. 

There are only a handful of objects whose cavities have been revealed in resolved images 
at multiple wavelengths covering a broad range. A special group of them are circumbinary 
disks (with the secondary being a stellar or sub-stellar object), such as GG Tau. The cavity 



in GG Tau has been found in the optical (Krist et al. 2005), NIR (Itoh et al. 2002), mm 



(Pietu et al. 2011), and CO line emission (Guilloteau et al. 1999), all at roughly the same 



position. The formation mechanism of cavities in circumbinary systems has been studied 
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by Artymowicz & Lubow (1994). It is thought that the gravitational interaction between 
the massive (sub-)stellar secondary and the disk naturally truncates the disk, and the cavity 
formed in this way exists in all disk components, including gas and dust at different sizes. 
The imaging at V band shown in Hashimoto et al. (2012) has put an upper limit ~ 30 — 50Mj 
for any possible companion at the radii of interest (assuming their age of the system and 
the (sub-)stellar object evolution model), which rules out the possibility of a stellar mass 
companion. Future observations with better sensitivity and contrast performance are needed 
to answer whether a sub-stellar companion exists (i.e. a brown dwarf). On the other hand, 



the possibility of gap opening due to multiple Jovian planets has been explored by Zhu et 



al. (2011); Dodson- Robinson & Salyk (2011). The general picture in that scenario, that a 



wide deep gap is opened while a small optically thick inner-most disk is left at the center, 
qualitatively agrees with the observations of PDS 70, though it is generally harder to directly 
image Jovian planet (s) in a bright protoplanetary disk. 



The dust filtration model has been proposed by Paardekooper & Mellema (2006) and 



Rice et al. (2006), in which the pressure maximum at the planet-induced gap outer edge acts 



like a filter, so that big grains are trapped but small grains penetrate into the inner disk. 
As a result, it is predicted that the cavity must be depleted of big dust grains, but with a 
significant amount of small dust particles left inside. PDS 70, with its heavily depleted cavity 
in small dust grains, is inconsistent with the dust filtration model alone. On the other hand, 



Zhu et al. (2012) proposed a dust filtration + dust growth model as a possible explanation for 



the transitional disk GM Aur. The grain growth and coagulation at the innermost disk may 
turn the small grains into big grains, resulting in a (pre-)transitional-disk-like NIR excess 
(see also Birnstiel et al. 2012). The dust and gas components of the disk are decoupled in 
their model and provide an explanation for both the moderate accretion rate of GM Aur and 
its strong near-IR deficit. Furthermore, this dust size dependent filtration model may explain 
the different gap properties between near-IR and sub-mm reported in Dong et al. (2012). 
However, unlike GM Aur, PDS 70 is a WTTS, suggesting a very low gas accretion rate. 
Together with the clear cavity revealed in the NIR imaging, a consistent picture emerges 
for PDS 70, that both the small dust and the gas are heavily depleted at the inner disk, 
and there is no clear evidence for decoupling between the two. On the other hand, without 
sub-mm/mm observations, we do not know if small and big dust are decoupled or not. 

Nevertheless, following the dust filtration + dust growth model, we can still provide an 



explanation for the PDS 70 class of objects. They may start as systems modeled in Dong et 



al. (2012), which do have a significant amount of small dust inside the cavity and the NIR 



cavity is not present. Later on, grain growth and coagulation, which happen at a faster rate 
in the inner disk due to high density and short dynamical timescale, gradually spread over the 
entire inner disk, eventually leading to the formation of a small dust cavity. If this scenario is 
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true, then these two kinds of (pre-)transitional disks are just at different evolutionary stages 
in their cavity clearing process — the ones in Dong et al. (2012) are at an early stage, while 
PDS 70 is at a later stage. Consequently, we should expect to see a statistical "time delay" 
between the two types of objects. This time delay is partly supported by the facts that 
PDS 70 is a WTTS, and that it is a relatively old system. This argument can be further 
tested if a large sample of objects in both categories with accurate age determination could 
be provided by future observations. Also, another possible strong evidence in favor of this 
scenario would be the detection of objects in the intermediate phase of this process, such 
as (1) objects with a "partially cleared" small dust cavity whose edge is between the center 
and the edge of the big dust cavity, or (2) objects with a radius-dependent cavity depletion 
factor for the small dust (probably smaller at inner disk), showing that the clearing process 
is moving outward. Furthermore, the byproduct of grain growth and coagulation — a slight 
enhancement of the big dust signal inside the cavity at later stage — may be observable 
using a sub-mm interferometer with very high sensitivity, such as ALMA (Dong et al.|2012). 



6. Summary 

We carry out a study of the disk+cavity structure for PDS 70, using radiative transfer 
modeling to fit both the observed SED and the SEEDS polarized scattered light image at H- 
band of this object. Good agreement with observations is achieved in our models. The disk 
has a giant cavity at its center with radius 65 AU. The small dust (sub-micron-sized) inside 
the cavity is depleted by a factor of ~ 1000, resulting in a low density but still vertically 
optically thick innermost disk, producing a pre-transitional-disk-like NIR excess. This heavy 
depletion is also needed to explain the surface density depression inferred from the scattered 
light image. The scale height of the small dust at the cavity edge is ~6 AU, constrained by 
both the SED and the image, which is consistent with the output mid-plane temperature of 
the disk. On the other hand, the total mass of the disk can be estimated only crudely to be 
on the order of 1O~ 4 M , due to the lack of sub-mm and mm data, and the degeneracy of 
the dust models for the big dust in the SED. This quantity, along with the cavity depletion 
factor for the big (mm-sized) dust, would be determined by future observations at longer 
wavelengths. 

Unlike most previously classified (pre-)transitional disk, the cavity in PDS 70 is identified 
in NIR scattered light, which only informs us of the depletion in the small dust but tells us 
little about the spatial distribution of the big dust. Pending the confirmation of the cavity in 
the big dust by radio interferometer observations, PDS 70 may be a prototype of its group, 
in which the cavity is seen in both dust populations. It is in clear contrast with the (pre- 
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)transitional disks discussed in Dong et al. (2012), where the small dust cavities were not seen 
outside the inner working angle of the scattered light images. The (pre-)transitional disks 
with or without a NIR cavity may be formed through different mechanisms (i.e. binary vs 
planets or grain growth), or they may just be caught at different evolutionary stage in their 
disk clearing process. Observational predictions for both mechanisms are made in Section [5j 
and more objects with multi-wavelengths observations in both categories are needed to reveal 
the nature of (pre-)transitional disks. 
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Fig. 1. — Opacity k for the dust models used in this study (showing the dust opacity only). 
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Fig. 2. — Surface density radial profile of the small dust S s , for models shown in Figure [5^,b. 
The horizontal dashed line indicates where the vertical optical depth is 1 at l/xm (character- 
istic wavelength of the stellar radiation). The fiducial model and model SCM1 have vertically 
optically thick innermost disks, and models SCM3 and SCM4 have optically thin innermost 
disks, while model SCM2 is on the margin. 
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Fig. 3. — Comparison of the polarized scattered light images at if -band between our fiducial 
model and SEEDS observation of PDS 70. Images show the surface brightness. The con- 
volved model image is the raw model image convolved by a HiCIAO PSF (Section [2]). The 
observed SEEDS image is not smoothed. The residual image is the convolved model image 
subtracted from the SEEDS image. All images are oriented in such a way that the far side 
of the disk (21° from west to north) is up. The mask at the center in the convolved model 
image and SEEDS image indicate a Qf!2 (radius) inner working angle. Labels are in unit of 
mJy/arcsec 2 (the residual image is on linear scale, while the other three are on log scale with 
the same color scheme). Our fiducial model matches the large scale characteristics in the 
observation well, although some local and asymmetric features are not reproduced perfectly. 
See Section [3] for details. 
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Fig. 4. — Comparison between observations and our fiducial model on the SED (a), SBRP 
of the scattered light image along major axes (b), and SBRP along minor axes (c). Ob- 
servational data is from paper I. In panel (a), the black dots are photometry data points 
(the vertical small bar indicates the error), and different components in the model SED are 
labeled. In panel (b) and (c), the black dots are the measured SEEDS SBRP, with error 
bars overplotted as vertical ticks. The individual observed SBRPs along the two directions 
of the semi-major axis are slightly different, and both are plotted as well. The fiducial model 
agrees well with the observations (particularly on the absolute scale of the surface brightness 
in the SEEDS image). In panel (c), due to the difficulty in measuring SBRP along the 
minor axes (Section [3J, the agreement is worse than in panel (b), but nevertheless the basic 
characteristic trend is still well matched, that the SBRP on the far side of the disk peaks at 
larger radii and decreases slower than on the near side, due to the back illumination. Model 
SBRP in (b) and (c) do not flatten out to a (constant) background noise at large radii, as 
in the observed SBRP. 
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Fig. 5. — Effects of the cavity depletion factor for the small (5 caVjS ) and big (5 caVj b) dust, on 
both the disk SED ((a) and (c)) and the scattered light image ((b) and (d), plotting SBRP 
along the major axes). The dotted line in panels (a) and (c) is the stellar spectrum and the 
black dots are photometric SED data point, as in Figure |4j The corresponding disk models 
are listed in the block "The Depletion Factor Inside the Cavity" in Table [TJ and described in 
detail in Section 4.1 The models in the SBRP plots ((b) and (d)) have the same line types as 
they have in the SED plots ((a) and (c)). The results show that as long as the innermost disk 
is vertically optically thick, the SED only weakly depends on 5 caVtS ; once <5 caVjS drops enough 
for the innermost disk to be optically thin, the NIR excess becomes sensitive to the amount 
of small dust there. The inner disk is brighter and the outer disk is dimmer in scattered light 
for a bigger <5 caVjS . On the other hand, both the SED and scattered light image are almost 
independent of 5 C av,b- 
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Fig. 6. — The same set of plots as in Figure |5j but for models BOM1 to BOM4-flat, showing 
the effect of total disk mass Md us t- The corresponding disk models are listed in the block 
"The Total Dust Mass of the Disk" in Table [TJ and described in detail in Section ^2 The 
dust model for the big dust in models BOM3-flat and BOM4-flat is the flat-big-dust, while 
the other models have the steep-big-dust model (Figure [T]). The total mass of the disk, which 
in PDS 70 concentrates in the big dust in the outer disk, needs to be better constrained by 
observations at long wavelengths, with which different big dust models can be distinguished. 
The photometry at the longest wavelength available (160/mi) agrees with a total disk mass 
~ 3 x 1O~ 5 M assuming the steep-big-dust model, or ~ 1.5 x 1CT'' 
dust-model. 
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Fig. 7. — The same set of plots as in Figure [5j but for models SOS1-4 and SCS1-4, showing 
the effects of scale height of the small dust, both outside the cavity (h°, (a) and (b) and 
inside the cavity (h°, (c) and (d)). The corresponding disk models are listed in the block 
"The Scale Height of the Small Dust" in Table [TJ and described in detail in Section 4.3 Both 
the IR excess at ~ 40/mi and the brightness of the disk at the cavity edge sensitively depend 
on the overall scale of h° (a higher cavity wall produces more MIR excess and scattered 
light), while if the overall scale is roughly the same, the radial dependence of h° does not 
have a big effect. Similar results are seen in K as well. 



